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Abstract 

We discuss a six dimensional mass generation for the neutrinos. Active neutrinos 
live on a three-brane and interact via a brane localized mass term with a bulk six- 
dimensional standard model singlet (sterile) Weyl fermion, the two dimensions being 
transverse to the three-brane. We derive the physical neutrino mass spectrum and 
show that the active neutrino mass and Kaluza-Klein masses have a logarithmic cutoff 
divergence related to the zero-size limit of the three-brane in the transverse space. 
This translates into a renormalisation group running of the neutrino masses above the 
Kaluza-Klein compactification scale coming from classical effects, without any new 
non-singlet particles in the spectrum. For compact radii in the eV-MeV range, relevant 
for neutrino physics, this scenario predicts running neutrino masses which could affect, 
in particular, neutrinoless double beta decay experiments. 
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1 Introduction and Conclusions 



Extra dimensional models have become quite popular in the past few years with the realisa- 
tion that they could exist at far lower scales compared to the Planck scale, Mpi, and thus they 
could have an impact at the TeV scale physics PQ 12 E] • Some of the significant implications 
could be for the gauge coupling unification jlj, solutions to the hierarchy problem [UlSj each 
having its own phenomenological impact like for example, collider searches jH] and precision 
measurements of standard electroweak parameters 

The existence of the extra spacetime dimensions at some scale closer or higher than the 
weak scale also gives much scope to build models related to fermion masses and mixings, par- 
ticularly for light neutrinos 0, qualitatively rather different from the traditional 4D seesaw 
mechanism |Hj. In this case, a singlet neutrino (right handed) could be allowed to propa- 
gate in the bulk whereas the SM (left-handed) neutrino would be confined to the 3-brane, 
leading to new possibilities of neutrino masses and mixing jHl HHl IE] • While this possibihty 
has been explored in detail both in model building as well as in terms of phenomenological 
analysis most of these studies have been confined to the case where there is only one 
additional dimension. It is necessary to extend these analysis for higher number of extra 
dimensions. This is because new features could arise in higher dimensional field theories 
which could give rise to different phenomenology. In particular, interesting 6D models have 
been constructed to break the electroweak symmetry through non-trivial Wilson lines [[13] , 
to guarantee the proton stability up to dimension fifteen operators J3] , to predict the num- 
ber of chiral generations to provide a dark matter candidate [12], to construct realistic 
GUT models ^7] etc. 6D models have also been constructed to reproduce the neutrino mass 
spectrum [TH]. 

In the present work, we will study the case of neutrino masses in a six dimensional 
model. In six dimensional models, it has been known for some time that orbifold compact- 
ification produces some peculiar properties such as 'tree level' renormalisation of coupling 
constants 120] whereas one-loop quantum effects in open string theories have a similar 
interpretation [21]. While this has been known for the scalar case for some time, we will 
study this property in detail for the case of neutrinos with a brane localized mass term on 
a T^/Z2 orbifold. We find that the physical neutrino mass (and also higher Kaluza-Klein 
modes) have a logarithmic divergence related to the brane thickness. We find this result by 
two different methods. First by diagonalizing the (infinite) mass matrix in the KK basis, 
along the lines of [3 dOl 122] • Second by considering the bulk propagation of the fields with 
appropriate boundary conditions due to the orbifold projection in presence of brane-localized 
operators [221 ■ 

We show that its interpretation is similar to the brane-localized scalar mass term discussed 
in and it can be renormalised in a similar manner by adding a neutrino Dirac mass brane 
counterterm. As a result, the neutrino mass runs, more precisely increases with energy, above 
the compactification scale -R~^, where for simplicity we consider the case of two equal radii. 
For radii in the eV to MeV range, the effect of this classical running has no counterpart in 
four dimensions, since it arises without the presence of any new particle charged under the 
Standard Model gauge group. This can be tested in processes with off-shell neutrinos, like 



1 



the neutrinoless double beta decay, where the increase in the neutrino mass at GeV energies 
enhances the amphtude of the process. 

The paper is organized as follows. Section 2 describes the example of a six-dimensional 
scalar field with four-dimensional (brane-localized) mass term in one of the orbifold fixed 
points, in the simplest orbifold compactification. We diagonalize explicitly the mass matrix 
and show that the physical masses (eigenvalues) have a logarithmic dependence on a UV 
cutoff, related to the thickness of the brane where the mass term was inserted. The same 
result is then obtained by solving bulk field equations with appropriate boundary conditions 
at the mass distribution, giving therefore a clear meaning of the ultraviolet cutoff in terms 
of the profile of the mass distribution. In Section 3 an off-shell analysis along the lines 
of [ini I2D1 suggests that the cutoff dependence has a natural interpretation in terms of 
renormalisation group running of the localized mass term, between the compactification 
scale and the ultraviolet cutoff. In Section 4 we pass to the case of interest in the present 
paper, four- dimensional localized active neutrino mixing via a localized Dirac mass term 
with a bulk six-dimensional singlet Weyl fermion. We show that the eigenvalue equation in 
this case is exactly the same as in the scalar field with localized mass discussed in Section 2. 
We then break the lepton number conservation by adding a brane-localized Majorana mass 
term for the bulk field and work out again the physical masses and their corresponding 
energy running. In Section 5 we discuss higher- dimensional neutrino oscillations and in 
particular the probability of regeneration of active neutrinos. We end up in Section 6 with a 
quick survey of the phenomenological consequences of this class of models coming from the 
fundamental scales and the running, which can be tested in particular in the neutrinoless 
double beta decay experiments. 

2 The Scalar Case 

We will discuss the 6D formalism in detail both in the scalar as well as the fermion case. We 
consider a T^/Z2 orbifolding in the two dimensional compact space. As a starting point, let 
us set the bulk mass to zero and add a localized mass term at the origin of the 2D compact 
space. The corresponding action reads'^ : 

S=^l d'xd'y ((9^^<|.)(9m$) - h $'5'(y)) , (1) 

where /12 > is a dimensionless coupling in the natural 6D units for which the scalar field 
$ has dimension two. The coupling /12 is localized at the origin of the compact space. The 
field equation is free in the bulk and has a delta function source at the origin 

9Af9^''$ + /i2$5'(y) = . (2) 

^We are using a (+,—,—,—,—,—) metric. The index M denotes bulk coordinates and runs from 
0,1,2,3,5,6, while fi — 0,1,2,3 denotes brane coordinates. We'll use either x^'^ or yi^2 to denote the 
two extra dimensions. Finally, 81^2 is a short-handed notation for dy-^ ^ 
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We will consider compactification on the orbifold T'^/'L2-, acting as the reflection {y 1,1/2) 
{—yi, —y2)- This orbifold has four fixed points as summarised below in Fig.(l). The fixed 
points and the corresponding Z2 transformations of the coordinates around those points are 
summarised as: 

yi -yi yi -Vi + 27ri?i yi ^ -y^ yi -yi + 2nRi 

1/2 -y2 2/2 -2/2 2/2 -2/2 + 27ri?2 2/2 ^ -2/2 + 27ri?2 (3) 

(0,0) Mi,0) (0,7ri?2) (7r/?i,7ri?2). 

In complex notation, the action of Z2 on the compact space is a two-dimensional tt rotation, 
^2(2/1 + ^2/2) = €.™{yi + iy2). Let us now proceed to study the KK spectrum. 




Figure 1: Z2 orbifold with a brane localized mass term for the scalar field at the origin. 



2.1 The KK approach 

If the scalar field $ is even under the orbifold action, it can be decomposed on a complete 
basis formed by the cosine functions: 

^{x,y) = YI (^i'2/2|fci,fc2)0(fci,fc2)(a;) (4) 

(fci,fc2)eJ 

with 

1 cos (^yi + ^^y2) 
(2/1, 2/2|fci, ^2) = p , , , , (5) 

The indices fci 2 belong to the set X 

X= {(O;0),(l...oo;0),(0;l...oo),(l...oo;l...oo),(l...oo;-oo...-l)} . (6) 
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The scalar action (P) then takes the following form after integration over the two extra 
dimensions 

(fci,fc2)e2: ^ 1 2/ \(fci,fc2)ex ^ / 

where 



is the naive (volume suppressed) four dimensional lightest scalar mass, which is typically of 
the order and slightly smaller than the compactification mass scale. The mass term of the 
4D action (0) is 

(fei,fc2),(pi,P2)e2: 

with the mass matrix given by 



2^2 / ^2 ^2 

The diagonalization of this mass matrix will define the KK mass eigenstates. 
2.2 Eigenvalue Analysis 

Let us now try to find the eigenvalues and eigenvectors of the mass matrix ()1U|] . The 
characteristic equation is given by 

M^^m = m'^m. (11) 

where represents the eigenvalues and ^ is the eigenvector in the basis {ki,k2)<^i 

de- 
fined in the previous section {'^{ki,k2) = {ki, fc2|^m) )• The matrix equation (fTT|) is equivalent 
to the infinite set of explicit equations for every (/ci, ^2) G X 

v^yn^^' = v/2V^ |^^2 _ ^ _ ^ ^^^^^^^^^ (^2) 
where we have defined 

^>' = V2 ^^vl>(,^ (13) 

(fci,fc2)ex ^ ^ 

The solution of the equations ()12|) is simply given by 



where is a normalization constant independent of {ki, ^2). Putting this solution back in 
the equation (fT^ for ki = k2 = 0, we obtain the eigenvalue equation we were searching for: 



^2 J2 rn^_ kl/Rl - kl/Rl ■ ^^^^ 

Note that, in this last equation, the sums have been conveniently rewritten from k^ = —00 
to +00 (see the identity (jH^ of the Appendix A). 

We want to find an estimate for the lightest solution, m^, of the eigenvalue equation. We 
are first going to evaluate the sums in ()15|) by an analytic continuation procedure in the mass 
eigenvalue m = ifi, justified by the absence of any poles in the complex m plane except for 
the Kaluza-Klein poles. Using the definition of the theta functions given in the Appendix, 
we can rewrite the eigenvalue equation as : 



2 - I dte->^'^9J^)eJ^, 

fOO 



dt e"^'* 



it \ „ ( it 



TlRfJ \7tK: 



?2 



-71R1R2 / — e Osi ^3 — ^ , (16) 



1/A2 t V ^ / V ^ 

where we have introduced an UV cut-off A noting that the integral is log divergent. In the 
following equations we consider for simplicity equal radii Ri = R2 = Rc, where Rc defines 
the compactification scale. We also rewrote, using properties of the ^3 function defined in 
the appendix, the last contribution in a more convenient way. 

Large Radii Limit: In the large radii limit, with respect to the cutoff A and by antici- 
pating that the lightest mass eigenvalue m < R~^ (we will discuss in section 3 the solutions 
corresponding to the massive KK excitations), the approximate solution of the above eq. lfT^ 
becomes 

= ^ ^ -T^Rl H^'Rl) + ■ (17) 

The second piece of the RHS of the above equation is the one corresponding to fcj = 0: it is 
the infrared contribution coming from the integration region i?^ < t < 00, whereas the first 
piece of the RHS, logarithmically divergent, is coming from the ultraviolet region of heavy 
states 1/A^ < t < R^. The first interpretation of the result is as follows. When the radii are 
small, we have the ordinary "volume suppressed mass" given by : 

-L.e = ^ • (18) 

On the other hand for large radii the logarithmic contribution from the UV has a significant 
effect and we get the renormalised contribution 
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This corresponds to the hghtest eigenvalue. There is also a tower of massive KK states, the 
corresponding masses and their running will be presented in Section 3 (see eq (HH)). 

Let us now move to understand the corresponding wave function. This wave function is 
easily obtained from the components of in the \ki, /c2)(fei,fc2)ex basis: 

^miy) = {yi,y2\'^m) = {yi^y2\h,k2){ki,k2\'^m) ■ (20) 

(fci,A:2)eX 

By using the eq.lfT^. this is evaluated to be: 

, \ _ ^_!)/^_ cos{kiyi/Ri + k2y2/R2) 

^"^^^^ " v/M^ ^ m2 - kf/Rj - kl/Rl ' ^ ^ 

fel,K2 = — OO 

where, again, we used the even properties of the cosine functions to transform the sums into 
sums from ki = — oo to +00. 

The same results can actually be obtained in a simpler way through appropriate boundary 
conditions, with special care in the fixed point where the localized mass lives. 



2.3 Understanding through Boundary Conditions 

As in 5D where we can reformulate orbifold projections in the presence of boundary operators 
in terms of boundary conditions for fields propagating in the fundamental domain we 
would like here to show that we can derive the eigenvalue equation by imposing suitable 
boundary conditions in 6D. 

In absence of any operator localized on the fixed points, the Z2 parity assignment for an 
even scalar field translates into the BCs at the fixed points 

= = . (22) 

In the presence of the localized coupling h2, the corresponding BCs can be read from the 
full equation of motion which, upon the KK decomposition $(x,y) = v^m(y) <Pm(yX), m 
being the 4D KK mass dfj,d^(f)rn{x) + rn^cprnix) = 0, reads 

To obtain the corresponding BCs at the origin, we first regularize the mass distribution by 
distributing it along a circle of radius e around the origin. The mass density is related to 
the localized dimensionless parameter h2 by 

/i2 = 27rem^. (24) 

With this linear mass distribution, the BCs along the surface yf + yl = {n is the normal 
to the surface) becomes 

dnip = -mlip (25) 
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and when e shrinks to zero, we obtain the new BCs at the origin 

dy2V = vKVi) ■ (26) 

We will now check that these BCs lead to the eigenvalue equation (fT^ derived earlier by 
the diagonalization of an infinite mass matrix. First, we note that the most general solution 
of the equation of motion ()2H|1 in the bulk takes the form 

^™(y)=Ar 2^ .2/02 2/02 — (27) 



m2 - kl/Rj - kl/Rl 

ki,k2=-oo 1/ i z/ z 



which verifies 



{df + di)ipm = -mVm +^f' Yl cos{kiyi/Ri + ^22/2/^^2 + «k) ■ (28) 

/Cl,fc2 = — 00 

In fl27j) . Ok are phases related to the position of the localised mass term. For a mass 
term localised at the origin, boundary conditions fix ctk = 0. Then, by using the iden- 
tity cos(%/i?) = 27iR S{y), we finally obtain 

i^f + ^|)ipm = -m^^m + ^7^^RlR2^f'5\y) . (29) 

The comparison of the boundary terms in (j^^ and (j^^ leads to the desired eigenvalue 
equation (fT3jl . If the mass distribution is regularised according to ()25|) . then the sum over 
the KK integers (fci, ^2) in ffT3j) is effectively cut for large momenta \ki\ < ki^^ax = -Rj/e- The 
size e is related to the UV cutoff defined in the previous section by A = 1/e. 

Alternatively, one can also derive the same eigenvalue equation by simply enforcing the 
modified BCs fl26|) . To this end, one needs to massage the expression of the wavefunction ()21|) . 
Actually, one of the sum can be performed by using the formula: 

cos(A;x) TT cosh a(7r — x) 
^ k"^ + a sinh air 



We obtain: 



ik2y2/R2 cosh ( {Ri/R2)^/k^ - m'^Rl {n-yi/Ri 

^(1/1,1/2) = -M'7iRiR2 H2 202 ^ 7 X 43« 

^ ^kl-vp?Rl sinh [{T^RxlR2)^kl - 



which leads to the expression: 

dy,^{y^ = 0, 2/2) = 2n^RiR2 M' 5{y2) . (32) 

Inserting this relation into the BCs at the origin, we are left again with the eigenvalue 
equation (fT3j) . The boundary condition in the second direction is satisfied in the same way. 
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3 Regularization Dependence and running 



The interpretation of the logarithmic divergence in the sum is actually rather clean, by 
realizing that the brane localized couplings, following ^Hj, do run in the sense of the four- 
dimensional renormalisation group. Our example is similar in spirit to the one analysed 
in [20]. It was shown there that the (euclidean) Dyson resummation of the scalar propagator, 
by including the brane localized mass insertions, in a mixed representation : 4d spacetime 
momentum space and 2d extra-dimensional coordinate space, is 

G'(g4,x,y) = L'(g4,x,y) - /i2 £'(g4,x, 0)^(54, 0,y) 

+hl Diq^, X, 0) D(g4, 0, 0)D(g4, 0, y) + • • • 

where D(g4,x, y) is the free scalar six-dimensional propagator of four-momentum and 
two-dimensional positions x, y. The propagator at the origin at the compact space sum in 
the compact space -D(g4, 0, 0) is defined via 

Jfe,0.0) ^ ^Dte.k) ^ 1 . (34) 

k fci,fc2 

where k = {ki,k2) are the Kaluza-Klein momenta in the transverse space. The pole of 
the resummed propagator (jH^ defines the physical four-dimensional masses qf = w?. One 
therefore gets, by going back in the Minkowski space 

47r^i?ii?2 _ ^ 1 

h2 ^ m2 - kl/Rl - kl/Rl ' ^ ^ 

which is nothing that the eigenvalue equation we derived earlier. The off-shell expression ()33|) . 
however, is more suitable for renormalisation purposes since by standard techniques it leads, 
for /i > 1/ Rc (in the case of equal radii, Ri = R2 = Rc, that we are considering here), to the 
renormalisation group flow of the /i2 coupling^ 

Indeed, an approximate evaluation of the propagator for different values of the four-dimensional 
momentum gives 
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V D(g4, k) ~ In — , for qj > 1/R, 
J^Diq^, k) ~ H^'R'c) + -^^^^^ ' for ql < l/Rl , (37) 



^Ref. [2Qj finds a /3 function l/(47r) for the orbifold [a = 1/2). The factor 2 of discrepancy with our result 
is coming from the different normahsation of the brane- locahzed couphng. We did normahse it on the torus, 
whereas Ref. (201 normahses it on the orbifold. The rescaling /12 — > /12/2 exactly accounts for the missing 
factor of 2 in the (3 function. 
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and therefore for > l/R"^ the propagator exhibits a logarithmic scale dependence typical 
of a renormalisation group running of the coupling /12. The coupling /i2 in the eigenvalue 
equation ()35|1 is interpreted as the value of the coupling at a high cutoff scale A, naturally 
of the order of the mass scale (higher- dimensional Planck mass or the string scale) of the 
microscopic theory. Following this interpretation, the eigenvalue equation (jSSl), solved under 
the physically relevant, for our case, situation m <^ ^/Ri (we assume Ri ~ R2 for simplicity), 
defines the pole mass. For energy scales ^ > I/-R12, the g| dependence can be renormalised 
away by an appropriate counter-term, which gives rise to the running of the mass eigenvalue: 

' 4.^ij.fljl + (38) 



m2(/i) \h2 2tt /i / h2{lJ,) 

The eigenvalue equation has other solutions corresponding to the massive Kaluza- 
Klein states. They can be easily computed in the leading approximation in /i2 ^ 1 by 
searching solutions of the form 

^kh) = {^ + ^](^ + ^h,h) , (39) 



" \Rl ^ Rl 
A straightforward computation gives 



Qj^ij is a degeneracy factor counting the number of solutions to the integer equation ll/Rl + 
l^/Rl = kf/Rf + A;|/_R|, with ki, ^2 = —00 ... -|- 00. For instance, in the case of equal radii, 
ci,i = 4 but Co, 10 = 12. 

Similarly to the lightest eigenvalue, the KK masses also run, but in this case the running 
occurs above their pole mass /i^ > m^^^ 

+ (41) 



4 The Fermionic Case 

Before deriving the mass formula, let us first describe our framework in more detail. We will 
assume that the Standard Model fermions are restricted to the 4D brane whereas a singlet 6D 
Weyl neutrino is free to propagate in the full six dimensional space^. The generalization to 
more than one singlet neutrino, needed in the three generation case, is straightforward. The 
six dimensional bulk action is simply 

S = j d^xd^y i^T^^VM^^ ■ (42) 

^In the presence of gravity or bulk gauge symmetries, anomaly cancellation in six dimensions gives 
nontrivial constraints on the 6D Weyl fermionic spectrum, which have to be taken into account in building 
a complete theory. 
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Notice that no Dirac nor Majorana bulk mass is allowed for a 6D Weyl fermion. 
The Z2 orbifold projection acts on ^ as 

Z2*(y) = ie-(-/2)r'r'^(-y) . (43) 

Decomposing the eight-component 6D Weyl spinors into two four-component 4D Weyl spinors 
(see Appendix B), the Z2 action becomes 

^ = ^ , Z2Ai(y) = Ai(-y), Z2A2(y) = -A2(-y) . (44) 

This allows us to write a brane locahzed couphng between the SM left-handed neutrino and 
the even fermion Ai: 

d^x {huLXi{y^O)H + h.c.) , (45) 



where H is the 4D Higgs field and h is the Yukawa coupling. The mass dimensions oi^ H 
are respectively 5/2,3/2 and 1. 

Written in two-component spinor notations, the 6D lagrangian takes the form : 

C = -iAi cr'' d^Xi - 1X2 8^X2 + Ai {d^ + ids) X2 

-X2 {d, - ide) Ai + g2{uLX, + ULXi)d\y) , (46) 

where we have now introduced 

92 = h{H) (47) 

the brane-localized 'Dirac' neutrino mass parameter which actually, in analogy with the 
scalar case, is a dimensionless parameter. We will study the phenomenology for two cases: 
(a) Dirac type neutrinos; (b) Majorana type neutrinos. Before proceeding further, we would 
like to point out that the equations of the motion are given by : 

ia^'d^Xi - {85 + ide) A2 = -92 '^iS'^iy) , 
^a^^^^X2 + {d^ + ide) Ai = , 
iai^d^UL = -92 Ai(y = 0) . (48) 

Using the identity on the Pauli matrices given in the Appendix B, one can easily combine 
these coupled first order differential equations to obtain an uncoupled second order differen- 
tial equation that is nothing but the Klein-Gordon equation in 6D: 

{d.d'^ - dl - dl) X,{x, y) + gl X,{x, y) ^^(y) = (49) 

and therefore our fermionic problem with brane localized Dirac mass term is reduced to the 
one of the bulk scalar field with brane localized mass term studied in the previous Sections. 
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4.1 Dirac Neutrinos 



Let us now study the mass matrix for the case of Dirac neutrinos. To this end we will 
follow the procedure first developed in the scalar case and we will expand the 6D fermion 
^ on a complete set of functions. Wc first write the 6D Weyl fermion in terms of a pair of 
two-component spinors, Ai^2- According to their Z2 parities, Ai^2 can be decomposed as: 



Ai(a;,y) = 



1 



y/2n'^RiR2 
1 



ECOs(/,-i(/i//?i + /,-2(/2//?2) , X 

(fci,fc2)ex ' 

^ sm{kiyi/Ri + k2y2/R2) C(fei,fe2)(^) • 



(50) 
(51) 



(fei,fe2)ex 

After integration over the two extra dimensions, the mass terms in the action become 

V2 



(fei,fe2)ex 



with 



m = , = . (53) 
V47r2i?ii?2 

We first notice that the phases in the KK masses can be absorbed by redefining the fields 

'C(fei,fc2)) while the phase in m can be removed by a redefinition of Ul. Thus, these phases do 
not have any physical effects and all the mass terms can be taken real. Next it is convenient 
to introduce the combinations (by convention, we define /l(o,o) = X(o,o) and -B(o,o) = 0) 



A 



(ki,k2) 



XikiM) + C(fei,fe2) 



V2 



and B 



(kiM) 



X(fcl,fc2) C(fcl,fc2) 

7i 



(54) 



The mass terms now write 



^{Pl,P2) 



(55) 



with 



M 



( ° 




m 


fh 






^ fh 








{P1,P2) 

fh 




^1 _|_ ^ 2 X 



(56) 



The mass matrix M. is real and symmetric. We will diagonalize is in a similar manner as we 
did for the scalar mass squared matrix. The characteristic equation defining the eigenvector 
associated with the eigenvalue m is 



= m 



(57) 
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This matrix equation is equivalent to the infinite set of equations 



m 

X 



rn{vL\'Em) + + -^(^{fci,fc2)|Sm) = "^(^(fci,fc2)|Sm) , (59) 



2 



^2 

m(z/L|Sm) - ■^{B{kiM)\-m) = m{Bi^ki,k2)\-m) ■ (60) 

Plugging back in eq. (j3H|) the expressions for {A(^f:i,k2)\'^rn) and (-B(fc^_fc2)|Sm,) obtained from 
eqs. (jSnHSni), we finally obtain the desired eigenvalue equation 

4.2 Majorana Case 

We now discuss the case of Majorana neutrinos. To the lagrangian discussed above (j46j) 
we want to add lepton number violating Majorana mass terms involving the Kaluza-Klein 
states of the bulk fermion \E'. Note that there are two ways to break the lepton number : 
(a) Breaking the lepton number on the brane by introducing on the brane a lepton number 
violating mass term for the singlet neutrino; (b) Breaking the lepton number in the bulk 
through a bulk Majorana mass for the singlet neutrino. We will directly present the mass 
spectrum of these two cases in the following. Let us first study the case (a) and let us add 
the lepton number violating mass term Mq on the brane 

j d^x d^y (Mo AiAi + h.c) 6\y) . (62) 

Notice that actually from a 6D perspective Mq has mass dimension —1, whereas after the 
KK expansion the physical mass parameter is Mq/{4tt'^RiR2). Similar to the Dirac mass 
considered previously, all phases in the Kaluza-Klein complex masses can be redefined away 
and have no physical meaning. By a straightforward generalisation of the previous diago- 
nalisation, we find the eigenvalue equation 



2 



m? + mMo/UTT^RM ^ w? - kVR^ - k^/R^ 

fcl,fe2 = — oo 



(63) 



By considering again for simplicity the case of two equal radii Ri = R2 = Rc and evaluating 
as before the double sum by keeping the leading IR and UV contributions, we find in the 
large radii limit 

—2 = ( 2^^ u + T- • (64) 

\g2 + mMo Air j 
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The natural interpretation is again in terms of the running of the physical mass 



(65) 



The case of the brane Majorana mass is the simplest but also the most problematic, since 
due to the double volume suppression the lepton number violation, is small. The case of 
the bulk Majorana mass M is more subtle. First of all, a standard Majorana, Lorentz 
invariant mass in 6D, ^'-''^ cannot be written for a Weyl fermion, since it mixes 6D Weyl 
fermions of opposite chiralities^. On the other hand, a Majorana mass term involving only 
a 6D Weyl fermion, of the form ^'^Tq'^ can be written, at the expense of breaking the 6D 
Lorentz symmetry which could be considered as being spontaneously generated by the vev 
of some vector field. In this case, it is not possible anymore to eliminate the phases in 
the KK masses by field redefinitions, even for real Majorana mass. Interestingly enough, the 
interplay between Kaluza-Klein masses and bulk Lorentz violating Majorana mass generates 
CP violation. This observation could be related to previous proposals to relate the CP 
symmetry to discrete subgroups of a higher-dimensional Lorentz group [24J. The eigenvalue 
equation in this case is 

2 °° -1 



X] ^ _ A/T _ /R, -L ^' ^^^^ 



47r^RiR2 ^ m-M - ikJRi + k2/R 

k\,k2=—oo 



5 Neutrino oscillations 

Generally, the gauge neutrino eigenstates z/j are related to the set of mass eigenstates z/, 
through a unitary mixing matrix U as 

vj = ^UfiVi, (67) 

i 

where the matrix U is extracted from the neutrino mass matrix. The probability of oscillation 
between two gauge eigenstates z/j and z/^ after a time t is given by 

Pf^At) = E \US'^S'^\' + 2 E (/7;,^;,,f/;^.f/;,,eW^^-^^W) , (68) 

i i>j 

where Ei = \/p'^ + mf is the energy of the mass eigenstate z/j. Of a particular interest for 
our case where the active neutrinos oscillates into bulk states is the survival probability after 
a time t, given by 

2 

(69) 

^E.D. thanks P. Holstein and S. Lavignac for discussfons on this issue. 
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We consider in the following the case of the Dirac neutrinos discussed in detail in Section 4. 
The matrix U in this case is (doubly) infinite dimensional and, from ()58|) - fl60|) its explicit 
form is 



(S»n|^l^(fcl,fc2)) = - 



1 



- ^kl/Rl + kl/R 
ffi 



{ki,k2) e J 



(A;i,A;2) e J\(0,0) 



mi + ^/kl/Rl + kj/R:^ 
The unitarity property of the U matrix determines the normalisation constant, Ni, to be 



(70) 
(71) 
(72) 



AT, 



ki,k2=- 
oo 

^1,^2 = — 



(m? - - kl/Rlf 

1 



[m. 



j - kl/Rl - kl/Rlf 



(73) 



where the eigenvalue equation (fT3j) has been used to arrive at the last expression. We notice 
that the last summation ()73|) is manifestly UV finite and therefore the only cutoff dependence 
arise implicitly through the running of the physical neutrino masses. The active neutrino 
survival probability, in the non-relativistic approximation, is then given by the formula ^ 



it) 



(fei,fe2)GX 



''(fcl.fc2)' 
2p 



fci,A,-2 = — 



2p 



[kiM) 



(74) 



The factor 2 in the first sum accounts for the fact that each massive state corresponds to 
a Dirac neutrino. An approximate estimation of the normalisation factors N(^ki,k2) gives the 
result 



(kiM) 



^2/ -^2 



(75) 



where rrii^kiM) P^^^ mass eigenvalue given by the eq. pUj) and Ck^,k2 are the degeneracy 

factors defined in Section 3. We can anticipate from (fTSj) some similarities and also some 
differences between the 5D and the present 6D case. Like in the 5D case [9j, for large KK 
masses Ni ~ ml and the active neutrino mostly oscillates into the lowest mass states. In the 
6D case, on the other hand, the degeneracy of the massive states is higher than in 5D and 
the decoupling of the massive states is slower than in 5D. 

Numerically too, we find similar behaviour. For a relatively small number of states, P{t) 
is always smaller than 1. In fact, even P(0) is smaller than 1, unless one sums over all 
the states. P{t) oscillates with t, though it is always remaining smaller than P(0). The 
oscillations become more rapid as we sum over more and more states. 

The brane Majorana case can simply be recovered from the previous expressions by the 
replacement fn^ + rriiMo. 



14 



6 Phenomenology 



3 




1 I , . , ^ , . , ^ , . , , , . ^ 

10-10 >^Q-e -10-2 10° 102 10"^ 

(GeV) 



Figure 2: Running of the neutrino squared masses with the energy. We have chosen A = 10^ 
GeV, g2 = I, and Ri = R2 = I eV~^. The running occurs between the cutoff scale and 
the compactification scale. Below the compactification scale, there is not running anymore. 

In this Section, we will detail some phenomenological consequences which could be of 
interest for neutrino masses and mixings. First of all, let us try to see the significance of the 
logarithmic running with some realistic numbers. Fig. |21 plots the running of the neutrino 
mass-squared with the energy. We have taken the radii to be of the order of 1 eV~^ with 
the fundamental Planck scale taken to be around 100 TeV, in order to recover the 4D 
Planck scale Mp through the standard relation RiR2Mf = Mp [T]. Note that these numbers 
are close to the present limits from cosmology on the size of two large extra dimensions [26j. 
The coupling g2 is set to one at the cutoff (microscopic) scale. We plot eq. (jHHj) and we 
see that the neutrino mass changes by at least a factor of 3 while running with the radius. 
Below R~^, there is no logarithmic running anymore and the masses keep their values at 
the compactification scale. This explains the horizontal line below R~^. We interpret this 
result as follows. The pole of the propagator as given in eq. (jH^ . defines the physical mass 
of the neutrino which is equivalent to the 'running' eigenvalue evaluated at R~^. As long as 
neutrinos are 'on-shell', they carry this mass. The 'running' neutrino masses on the other 
hand varies with the off-shell momentum. Experimentally, the effect of this 'running' could 
be possibly seen in processes where neutrinos are off-shell. The best probed process in this 
regard is neutrinoless double beta decay which involves lepton number violating neutrino 
masses. We note only that as long as the lepton number is violated by Majorana masses 
and the running of the Dirac masses does occur, this is enough to have observable effects 
due to the running between the pole mass and the off-shell momentum, of the order the 
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nucleon mass, in the neutrinoless double beta decay diagrams. A detailed analysis of this 
case is beyond our goal here. On the other hand, irrespective of whether Icpton number 
is violated or not, this running should always be considered while constructing models of 
fermion masses in 6D. This is especially true in models where large radii are considered. 
One can interpret that the 'bare' couplings in the lagrangian given at the scale A would get 
modified as they 'run' towards their physical values. While most of the present discussion 
has been concentrated for the case of a single generation, we can easily extend the above 
analysis for more than one flavour. First of all, it is useful to recall that in order to generate 
three non- vanishing masses for the three active neutrinos, we need to introduce three distinct 
bulk sterile Weyl fermions a-nd therefore the brane-bulk mixing mass term becomes a 3 x 3 
matrix. The obvious change in the RGE evolution of the dimensionless coupling /i2 = 92 is 

= ^ ih2)ac{h2U , (76) 

whereas the 3x3 lightest eigenvalue mass matrix has the RGE evolution 

i^r^) = 47r2i?ii?2 [(1) +^5^^in-]^A7r^R,Rj-^) . (77) 

\m^{U,)/ab \ho/ab ZTT it \ho{ll)/ab 



The classical running then modifies the fiavour mixing present at the scale A by renormalising 
it in evolving towards low energy. Since in order to have sizable running the active neutrinos 
are required not to be much lighter than the sterile neutrinos, the oscillations into the bulk 
sterile states are constrained by the existing experimental data. The two flavour active- sterile 
oscillations can be suppressed by choosing a sufficiently small coupling /i2. For example, for 
/i2 = 0.05 we find (for A = 100 GeV and summing up to first 200 states^), that the survival 
probability of the active neutrinos can stay around 70%. The survival probability obviously 
decreases as /12 increases. When other active fiavours are added, since mass differences 
for active neutrinos Am^^^^, Am^^^^^ <^ ~ eV"^, for small values of the coupling (72 
and for all KK states the mass differences m^. — ''tt-'kr ^ ^^^ej^m' ^"^^ej^t- ^^^^ case, the 
oscillations into sterile neutrinos are subdominant with respect to active neutrino oscillations. 
Obviously, the above results are qualitative and a full phenomenological analysis is required. 
We postpone it for a later work. 

It is more useful to write the running neutrino mass directly as a function of the physical 
pole mass m^. In the case of Dirac neutrinos and for one generation, for simplicity, we get 

9. ^ 9 l + #ln(Ai?e) , , 

m'M = rnl ^V^-r^ • (78) 

1 + iM^) 

This running is a small effect for small couplings g2- For example, for radii ~ eV and in 
the normal hierarchical neutrinos scenario, the coupling g2 for the electron neutrino has to be 

^For these values, the running of the neutrino masses from the cut-off scale to eV could be as 

large as 16 %. 
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of order 10~^. In this case, the running ()78|) in the energy range eV and GeV is too small to 
have observable effects. In the inverse hierarchy (or almost degenerate) neutrinos scenario, 
however, for g2 ~ C>{1), the effect is large, as already stressed. This fits well with the fact 
that neutrinoless double beta decay amplitude, proportional to the effective Majorana mass 

rriee = Yl ^ea^^ . (79) 
a 

where a are the three active neutrino states, increases with the absolute value of the neutrino 
masses. 

The question then remains about other physical effects of the sterile neutrinos. For 
the large radii, Rc ~ 0{eV^^) — 0(keV^^), we should expect a large number of very light 
sterile states. The effects of the virtual sterile states in loop processes involving neutrinos, 
in astrophysics and cosmology need a separate careful analysis, which is very important for 
the validity of the models we proposed in this paper. This is beyond the goals of the present 
work. We notice here, however, that the running of the neutrino masses is valid also for 
smaller values of the compact radii, even if in this case the effect of the running becomes 
much smaller and it will be harder to detect experimentally. 



Appendix A: Jacobi function 

The Jacobi function 6*3 used in the text is defined as 

e3{z\T)= J2 e^^'^'^'^^^^ . (80) 



fe=— 00 

A useful modular transformation property is 



^3 (^'v) =v^e-^'^' ^3(^,r) 



^1) 



The asymptotic limits of ^3 that we used are 



hm d^{z\T = it) = l , lim 6'3(z|t = it) ~ ^ e""^ . (82) 



t— +00 



We denote ^3(0|r) by e^{T). 

We can rewrite the eigenvalue equation (|THjl by using the Jacobi 6^ function. Indeed, by 
introducing the Schwinger proper time parameter t through 



1 

I 



00 

tA 



dt e"*^ (83) 
and using the analytic continuation prescription m = ifi, we find 



ki,k2=—oo 
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Let us also mention the identity 

2 °° 
E ^4:^/(^1'^^)= E /(^i'^^)' (85) 

(fci,fc2)GX fcl,fc2=— oo 

valid for any function / such that f{—ki, — ^2) = /(^i, ^2) and where X is the set defined 
in eq.®. 

Appendix B: 6D Dirac matrices 

For completeness, we give in this appendix the convention about spinors and Dirac ma- 
trices used throughout the paper. We have mainly followed the conventions of Wess and 
Bagger 

In 6D, the Dirac matrices are 8x8 and they can be easily constructed from the Dirac 
matrices in 4D: 

°^ r-f ° 't) a„dr-f5 7^ (86) 



with 

7''=(-°. and ,= = n»...7'=(7° | (87) 

where a'^ and a'^ are the usual Pauli matrices 

a' = -I2, = 

= a", a' = -a\ = -a\ = -a\ (89) 



01\ 2 o I 

1 = U = -1 



A famous relation about the Pauli matrices that is useful to link the fermionic equation of 
motion to the scalar equation of motion is 

a^CT'^ + (T^a^ = 27]'^'' and a^a^ + ctV^ = 27^'^^ (90) 

The chirality matrix in 6D is simply defined by 

r^ = r°...rW=(;' _° , ) . (91) 

Therefore a 6D Weyl fermion is written in terms of a pair of four-component Weyl fermions 

^={\[], (92) 



2 



with 

m = -r^^ ^ Ai = -7^Ai, A2 = 7^A2 . (93) 
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The four-component Weyl fermions can be themselves conveniently written as two-components 
spinors (we'll use the same notation for a four and two component spinor, the distinction 
should be clear from the context) 

\ , f \2a 



Ai = ( 1 , A2 = ( ) . (94) 



The dotted and undotted indices of a two-component spinor are raised and lower with 



the 2x2 antisymmetric tensors e^/j = io"^^ and e^^ = icr^A and their inverse e"^ = ~^<^a/3) 



ap 

^2 = ^"^>^2p and Ai^ = e.^Af . (95) 
Note also the adjoint relation: 

{\lr = K ■ (96) 

Finally, A1A2 and A1A2 denote the two Lorentz invariant scalars: 

A1A2 = A? Asa and A1A2 = XiaK ■ (97) 
These products are symmetric 

A1A2 = A2A1 and A1A2 = A2A1 . (98) 
The charge conjugation matrix in six dimensions, satisfying 

^pM^-i ^ -{V^f , (99) 

can be expressed in terms of the 4D one as 

C ^ ( ) . (100) 

where 

C4 1^ Cl^ = -(7'')^ ■ (101) 
In the chiral representation of the Dirac matrices that we are using, 

C.-7V^(*f _t,) (102) 

It can be easily checked that 

Cr^ = -r^C and & ^ C . (103) 
The charged conjugated spinor W is defined through 

= CW . (104) 
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It can be easily checked that for the 6D Weyl spinor satisfying eqs. ()93|). we have (in four 
and two component notations): 



= , = , 

^T^* = t (A*C4A2 + X\CAi) = - (AiAi + A2A2) 
^cp6^ = -A^C4A2 + A*iC4Ai = -AiAi + A2A2 . 



(105) 
(106) 
(107) 
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